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ABSTRACT 


The  fundamental  mechanics  of  deiaminution  in  fiber  composite  laminates  is 
studied.  Mathematical  formulation  of  the  problem  is  based  on  recently  developed 
laminate  anisotropic  elasticity  theory  and  interlaminar  fracture  mechanics  con¬ 
cepts.  Stress  singularities  and  complete  solution  structures  associated  with 
general  composite  delaminations  are  determined.  For  a  fully  open  delamination  with 
traction-free  surfaces,  oscillatory  stress  singularities  always  appear,  leading  to 
physically  inadmissible  field  solutions.  A  refined  model  is  Introduced  by  con¬ 
sidering  a  partially  closed  delamination  with  crack  surfaces  in  finite-length 
contact.  Stress  singularities  associated  with  a  partially  closed  delamination 
having  frictional  crack-surface  contact  are  determined,  and  are  found  to  be  dif¬ 
ferent  from  the  inverse  square-root  one  of  the  frictionleas -contact  case.  In  the 
case  of  a  delamination  with  very  small  area  of  crack  closure,  a  simplified  model 
having  a  square-root  stress  singularity  is  employed  by  taking  the  limit  of  the 
partially  closed  delamination.  The  possible  presence  of  logarithmic-type  stress 
singularity  is  examined;  no  logarithmic  singularity  of  any  kind  is  found  in  the 
composite  delamination  problem.  Numerical  examples  of  dominant  stress  singular¬ 
ities  are  shown  for  delaminations  having  crack-tip  closure  with  different  fric¬ 
tional  coefficients  between  general  8^  and  02  graphite-epoxy  composites. 
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i.  xmosocrxoB 


Oilulutioo  hu  bm  a  problem  of  significant  con  com  la  cho  rtllablo 
design  and  analysis  of  advanced  fiber  composite  Laminates.  Separation  of 
composite  laminae,  caused  by  high  local  Interlaminar  stress  aad  low  strength 
aloeg  the  ply  Interface,  can  result  In  destruction  of  load  transfer,  reduction 
of  stiffness,  and  loss  of  •true' >'ral  integrity,  leading  to  final  structural 
and  functional  failure.  From  the  mechanics  point  of  view,  de lamination 
Involves  Initiation  and  growth  of  macroscopic  cracks  between  dlsslmiler, 
strongly  anisotropic  solids.  A  rigorous  mathematical  study  of  delaad nation  Is 
recognised  to  be  difficult,  especially  In  a  finite -dimensional  fiber  composite 
laminate.  The  complexities  Include  the  inherent  crack-tip  singularity,  the 
effect  of  anisotropy  of  each  constituent  fiber  lamina,  and  the  abrupt  change 
of  stiffness  or  ply  orientation  through  the  laminate  thickness  direction*  la 
addition,  the  three-dimensional  state  of  stress  aad  deformation  associated 
with  the  coeposlte  delamination  always  gives  rise  to  a  combined  opening  (node 
1),  la-plane  shearing  (mode  11),  and  out-of-plane  tearing  node  (node  III) 
fracture,  which  render  the  problem  mathematically  intractable  in  many 
cases.  The  mechanics  of  de  lamination  in  fiber  composite  laminates  la, 
therefore,  not  only  of  significant  academic  interest  but  of  practical 
importance.  In  this  paper,  the  first  of  two  articles  la  a  row,  the 
fundamental  nature  of  stress  singularities  end  associated  field  solutions  for 
a  da  land,  nation  la  a  fiber  composite  laminate  are  Investigated. 

(bring  to  the  aforementioned  complexities,  studies  on  am  interface  crack 
between  dissimilar  anisotropic  materials  have  been  limited.  Gotoh  [1)  appears 
to  be  the  first  to  examine  the  two-dimensional  problem  of  partial  debonding 
between  dissimilar  anisotropic  plates  under  a  plane  stress  condition. 

Clements  (2]  has  used  Seroh's  approach  (3)  to  study  the  problem  of  an 
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interface  crack  between  two  generally  aalaotroplc  half -spaces.  Wlllia  (4]  haa 
also  conducted  a  two-diaensional  sc  rasa  analysis  of  a  crack  on  cha  plana 
lncarfaca  of  two  bonded  dlssiailar  half-spaces.  Tha  analysis  haa  bean 
coabined  with  cha  usual  local  fora  of  Griffith's  virtual  work  arguasnt  to  giva 
a  failure  criterion,  involving  a  stress  concentration  vector  and  specific 
surface  energy  of  tha  bonded  interface*  All  of  tha  asynptotie  solutions 
obtained  in  11,2,4)  have  an  oeclllacory  displacsnent  field  that  notarial 
lnterpenecracloo  on  either  side  of  the  crack  surface  la  predicted.  Slallar  to 
thoee  found  in  the  solutions  for  sa  lntsrface  crack  between  dlssiailar 
isotropic  asterlals  1 3-9),  these  physically  unreasonable  rssults  have  led  to 
the  arguasnt  of  solution  lnadalsslblllty  for  the  crack  problea  in  dlssiailar 
anisotropic  aedla.  Ito  correct  the  unsatisfactory  faature  of  oscillatory 
stress  singularity.  Hang  and  Choi  (10)  have  recently  reconsidered  the  problea 
of  an  Interface  crack  between  dlssiailar,  strongly  anisotropic  flber-coaposlte 
half-spaces  by  introducing  a  partially  closed  Interface  crack  nodal,  in  which 
tha  crack  is  not  coapletely  open  and  that  its  surfaces  are  in  frictionless 
contact  near  the  tip.  The  foraulatlon  leads  to  a  singular  Integral  equation, 
which  is  solved  nuaarlcally.  Nuasrleal  results  froa  this  refined  aodel  (10) 
exhibit  an  Inverse  square-root  stress  singularity  and,  therefore,  physically 
asanlngful  fracture  aschanies  paraaaters  can  be  defined  consistently  with 
those  In  fracture  probleas  of  homogeneous  aatsrlals  (11-13)  sad  in  the  aodel 
given  by  Coanlaou  (14,13)  for  an  interface  crack  between  two  isotropic 
asdla.  Moreover,  significantly  global  crack  closure  has  been  found  (16)  for 
an  interface  crack  between  dlssiailar  anisotropic  elastic  half-spaces 
subjected  to  aixad-aode  loading — a  situation  that  is  generally  experienced  by 
a  de lamination  in  finite  dimensional  fiber  composite  laainatas. 
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in  this  paper,  we  eaploy  Lekhnltskli's  eoeplex-varlabla  acraaa  potentials 
(17)  la  conjunction  with  «  eigenfunction  expansion  aathod  to  tualnt  tha 
Mchsnlca  sad  tha  aathaaatlcal  tolutloa  structure  for  a  dalaalaatloa  with 
frictional  crack-tip  closure  la  a  composite  laalaate.  Based  on  tha  general 
solutloa  structure  dateralasd,  an  advaacsd  auasrlcal  aathod  uslog  singular 
finite  aleasnts  la  than  developed  to  study  tha  detailed  dalaalaatloa  behavior 
la  finite  dlasnsloaal  fiber  coaposlte  laalnates  with  any  arbitrary 
coablnatlons  of  laalaatloa  paraasters,  geoastrlc  variables,  sad  crack 
dlasnslons.  Owing  to  space  Hal  tat  ion,  the  auasrlcal  aathod  sod  the  detailed 
coaposlte  delaal nation  behavior  are  reported  la  an  accoapaaylag  article  (18). 

In  the  next  section,  tha  problea  definition  sad  basic  ass uapt ions  era 
stated.  Basic  laalaate  aalsotrople  elasticity  equations  and  foraulatloa  of 
the  coaposlte  dalaalaatloa  problaa  are  Introduced  la  Section  3.  General 
solution  structures  for  asyeptotic  stress  sad  displacaaant  fields  are 
obtained.  Stress  singularities  associated  with  an  open  lnterlaalnar  crack  and 
with  a  partially  closed  delaal  nation  tip  with  frictional  crack-eurfaca  contact 
are  determined  respectively  In  Section  4.  Influences  of  frictional 
coefficients  on  dalaalaatloa  stress  singularities  are  exaalaad.  A  slapllfled 
aodel  for  a  delaal  nation  with  a  very  saall  area  of  crack-tip  cloeure  Is  also 
Introduced.  The  possibility  of  existence  of  additional  singularities  in 
logarithalc  fores  In  hoeogeneous  sad  particular  solutions  Is  Investigated. 
Results  are  presented  for  delsal nations  with  different  local  crack-surface 
traction  boundary  conditions  In  coaposlte  laalnates  containing  various  fiber 
orientations.  The  eigenvalues  and  associated  stress  singularities  obtalnedi  In 
this  study  provide  the  aost  fundaasntal  lnforaatlon  on  coaplete  solution 
structures  of  delaal  nation  stress  and  deformation  fields,  sad  establish  a 
basis  for  foraulatloa  of  the  singular  finite  elements  used  In  the  next  paper 
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(18)  Co  a  tody  cha  datailad  dalasi  oat loo  behavior  la  flalta  dlaaaaloaal 
coapoaltaa  with  gaaaral  laalaacloa  aad  gaoaatrlc  variables. 
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Tha  problem  conal  -  ;  -j  here  1*  i  eeapoitu  Ualuu  (fig.  1)  compose d  of 

unidirectional  fiber-reinforced  pllaa  of  uniform  eblekaaaaaa.  h|,  hj ,  hj . 

hjj.  Tha  eoapoelte  baa  a  flalca  diaaaaloa  with  a  width  equal  to  2b.  for 
siapliclty  and  without  loaa  of  generality,  wa  raatrlet  our  attention  to  tha 
eaaaa  of  symmetric  eoapoalta  laainataa  with  fiber  orleatatlooa 
( e&/e2/03/ • • • •  nf  thicknesses  are  alao  symmetric  with  reepect  to 
tha  *-*  plana ,  l.a.,  for  each  ply  abowa  tha  xi  plena  (y  >  0).  there  exists  an 
Identical  ply  with  tha  ease  ply  thiekaaaa  below  tha  w  plena  (y  <  0). 

Oelaai nation  with  a  length  a  la  eaauaad  to  occur  in  the  fora  of  an  Interface 
crack  between  dlaeialler,  strongly  anisotropic  flbar-ralnf oread  eoapoalta 
laaJnae  with  fiber  oriantatlona  •a  W 

The  eoapoalta  laalnata  la  eaauaad  to  be  subjected  to  tractions  acting  in 
planes  normal  to  tha  s-axls  and  distributed  uniformly  along  the  s-axis  without 
variation.  In  tha  case  that  tha  finite  dimensional  composite  laalnata  baa  a 
finite  length,  axial  Loads  and  moments  are  eeeuaad  to  act  on  the  ends  of  the 
eoapoalta  body,  the  eoapoalta  laalnata  la  further  assumed  to  be  sufficiently 
long  that  in  tha  region  away  from  tha  ends,  and  affects  are  negligible  by 
virtue  of  tha  Salat  tenant  principle.  Consequently,  tha  components  of 
stresses  in  tha  laadnata  are  Independent  of  tha  s-exls.  Tha  special  case  in 
which  all  components  of  stresses  and  displacements  la  tha  eoapoalta  are 
Independent  of  tha  a  is  well-known  as  the  generalised  plana  deformation 
problea  (17). 

The  objectives  of  this  study  are  to:  (1)  establish  a  mathematical  basis 
for  tha  aachanlca  of  da  lamination  baaed  on  leal  note  elasticity  theories  and 
interlaminar  fracture  mechanics  concepts;  (2)  determine  stress  singularities 
and  associated  eolation  structures  for  eoapoalta  de leal nations  with  different 
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local  crack -tip  da formation  cool lgura t loaa ;  (3)  obtain  aoyaptotlc  acraaa  aad 
doforaacloo  f lalda  governing  cha  fuadaaantal  behavior  of  da laal  nation;  and 
(4)  study  cha  Inf luoncaa  of  varlooa  loaf  not  loo  aad  notarial  variation  ouch  aa 
flbor  or ton tat loo  aad  crack -ourfaca  frlctloaal  caofflelooca  oa  cha 
do laal nation  ocrooo  singularities. 


3.  imma*  luincm  auu  rot 


Dm  dtwlopatot  of  the  mechanics  of  coapoelte  delamloacloa  Is  based  on 
recently  established  choorloo  of  anisotropic  lulaatt  elasticity  (19,20,21) 

•ad  tract an  mechanics  concepts  of  Interface  cracks  between  dissimilar , 
strongly  anisotropic  composites  (4,10,16,21).  la  this  sac t loo,  governing 
partial  dlffaraatlal  equations  for  coapoelte  laminate  a las tl city  problcas  are 
established  first,  baaed  oa  Ukhaltskll's  comp lax-var labia  stress  potsatlals 
(17).  General  eolations  for  the  leal  oats  elasticity  problea  with  lotarlaaloar 
cracks  are  latrodocsd.  Stress  singularities  associated  with  a  coaposita 
dslaalaatloo  hawing  homogeneous  local  boundary  conditions  are  defined. 

Solution  structures  of  asynpeotlc  stress  sod  dlsplacsnsnt  fields  are 
constructed  for  a  dslaalaatloo  between  dlsslallar  geoaral  fiber  coaposita 
lanlnae.  Additional  terns  of  logarithmic  forms  la  the  homogeneous  and 
particular  solutions  for  the  coaposita  da leal  oat loo  problem  era  examined. 


3. 1  basic  agnations 

tbs  fundamental  mechanics  of  delamloatlon  la  a  flber-relnforced  coaposita 
laminate  any  be  studied  from  the  schematics  illustrated  la  Fig.  1.  The 
constitutive  equations  of  each  flber-relnforced  composite  lamina  with 
recti  Honor  anisotropy  of  e  general  form  la  the  structural  (x-y-x)  coordinates 
are  denoted  by  pane  rail  sad  hooka's  law  la  contracted  ootetloo  as 

Cj  ■  Ijj  dj  (1,J  •  1 , 2 , 3, . . • ,6) ,  (1) 

where  the  repeated  subscript  indicates  emanation,  and  S^j  is  a  coapl lance 
tensor,  the  engineering  strains  t^  la  lq.  1  are  defined  la  a  Cartesian 


coordinate  system  by 
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3u 

el  ■  ex  "  IP 


e,  ■  e  -  ii, 
2  y  W 


3w 

e3 "  S  ■  yjT’ 


3w  3v  3w  3u  „ _ _  3v  ,  3u  ,,, 

e4“V’a7+3T’  e5  ■  Yzx  "  37  +  aP  e6  "  V  “37 +  37*  (2) 

where  u,  v,  and  w  are  dlsplacemenC  components.  The  stresses  are  defined  in 
an  analogous  manner  in  the  Cartesian  coordinate  system.  For  a  composite 
laminate  in  the  aforementioned  loading  condition,  mathematical  formulation  for 
this  class  of  elastostatic  problems  can  be  made  using  the  well-known 
Lekhnitskii  complex-variable  stress  potentials  f j7],  F(x,y)  and  7(x,y), 


defined  as 


x  3y2 


32F 

a  - - , 

y  Sx2 


xy  3x  3y* 


3?  of 

Tyz  "  3x*  r cz  3y* 

Following  the  same  procedure  in  [17,19],  u«  can  easily  obtain  the  following 
system  of  partial  differential  equations  for  each  anisotropic  composite 


lamina : 


L3F  +  4T  "  “  **4  +  A1S34  "  A2S35’ 


h4F  +  Ljf  -  0, 


where  L2,  L3  and  L4  are  differential  operators  of  the  second,  third,  and 

fourth  orders  which  have  the  form: 

L  «  S  !L  -  2S  32  +  S  *L 

2  44  3x2  2S«  3x  3y  +  S55  ,y2 ' 


L3  "  "  *24  T7*  (*25  +  *46*  'V"-  *  (*14  +  *56*  ' *  *15  !  3 
3x3  3xz  3y  3x  3yz  3y3 


3x  3y2 


L4  "  S22  —  -  2S26 


3x3  3y 


+  (2S12  +  S66 


2 - 2S  — - - +  S  ,  (5) 

3x2  3y2  3x  3y3  ty1* 
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t-  .  . 

Or  FOOii  ^  i  / 


^lj  ’  Sij  "  Si3  S3j/S33 


(i,j  -  1,2, 4, 5, 6). 


The  constancs  Aj  and  A2  In  Eqs.  4(a)  and  4(b)  characterize  the  bending  of  the 
composite  body  in  the  x-z  and  y-z  planes  respectively,  and  A4  is  the  relative 
angle  of  rotation  about  the  z-axis. 

Assume  that  external  tractions  on  the  lateral  surface  of  the  composite 

cross  section  are  given  as  T  ,  T  and  T  ,  The  boundary  conditions  on  the 

x  y  * 

lateral  surface  9B  are  as  follows: 

a  n  +  t  n  -  T  , 
xx  xy  y  x* 

Txy  nx  +  ay  ny  "  V  (7) 

Txz  nx  +  Tyz  ay  "  V 

where  n^  are  directional  cosines  of  the  bounding  surface  9B.  The  conditions 
at  the  ends  of  the  composite  have  the  form: 

//  *  dx  dy  -  //  x  dx  dy  -  0,  //  a  dx  dy  -  P  , 

B  “  B  yZ  B  *  Z 

//  y  dx  dy  -  Mj,  //  az  x  dx  dy  -  Mj,  (8) 

5  B 

//  <TyZ  X  "  XXZ  y)dX  ^  "  Ht* 

B 

where  the  Integrals  are  taken  over  the  entire  area  of  the  cross  section  B,  and 
Pz,  Mp  Mj,  and  Mj  are  applied  force,  bending  moments,  and  twisting  moment  at 
ends  of  the  composite,  respectively. 


3.2  General  Solutions 

The  general  solutions  for  the  governing  differential  equations  have  been 


shown  (17]  to  have  the  fora  as 
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rO'-' 


6 

f  -  I,  W  *  Po-  <» 

k-1 


6 

"  ■  J,  \  Fk‘V  +  V  <“> 

k-1 

where  the  coeplex  variables  Z^  are  defined  asZ|l-x+tijty;Fo  and  are 
particular  solutions  of  the  nonhomogeneous  system;  the  prime  (')  denotes 
differentiation  of  the  analytical  functions  F^  with  respect  to  their 
arguments,  and  Mfc  ere  roots  of  the  following  algebraic  characteristic 


equation: 

VH)12(H)  ’  *3<,l)  “  °*  UU 

with  l2(p)  -  S55p2  -  2S45p  +  S44,  (12a) 

l3(p)  -  S15p3  -  (Su  +  S56)p2  +  (§25  +  ?>46)p  -  ^24>  (12b) 

14(P)  -  V*  "  2V'  +  (2312  +  366^2  -  2V  '  V  (12c> 

The  in  Eq.  10  are  complex  numbers  equal  to 

\  *  "  *3*,lk^*2*Mk*  "  "  ty^)/^^).  03) 

We  now  choose  the  form  of  ^(Z^)  as 

Fk(Zk}  “  °k  z£*2/l<6+l)<S+2)l,  (14) 


where  C^.  and  &  are  arbitrary  complex  constants  to  be  determined.  Substituting 
®q.  14  Into  Eqs.  3,  9  and  10,  we  obtain  the  homogeneous  solutions  for  stress 
and  displacement  components  in  polar  coordinates  (r-$-z)  as  follows: 
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k-1 


“rh>  ■  +  Vs  V 


.5+1 


7T  V5+1  , 


k-1 


“i1*’  ■  .  S.lCk“7kC/(W'>  +  Vs  Sc  ^+1«5+l»-  »6> 


k-1 


.(h) 


♦  Vs  V 


rr  sr5+l , 


k-1 


where 


Zk  -  r(el*  +  +  Xk), 


Xk  -  (1  +  iuk)/(l  -  iMjj). 


(17a) 

(17b) 


The  coefficients  «lk(*)  (1  -  1,2,3 . 8;  k  -  1,2,3)  are  known  functions 

of  ♦,  n^,  uk,  and  ,  defined  In  Appendix  1. 

The  complete  laminate  elasticity  solutions  for  the  composite  mechanics 
problem  can  be  written  as 

°i  “  tflh)  +  °1P)*  (i  "  ^•M.S.b),  (18a) 

UJ  "  UJh>  +  UJP>’  <J  “  1,2,3)’  (18b) 

where  o^p^  and  ujp^  are  particular  solutions  associated  with  the  loading 
condition  of  each  Individual  case  studied.  The  expressions  for  and  o£p^ 


can  be  obtained  as 
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r  k 

«<h>  “  '  S3j  0,jh>/S33*  (l8c) 

end 

o<p)  -  (A^x  +  A2y  +  A^  -  S3j  o<p)/S33  (j  -  1,2, 4, 5, 6).  (18d) 

3.3  Asymptotic  Stress  sad  Displacement  fields 

Using  Eqs.  15  and  16  and  applying  local  homogeneous  traction  boundary 
conditions  on  crack  surfaces  3BC  (Fig.  2)  and  interface  continuity  (notching) 
conditions  along  3Bj,  we  obtain  a  syatea  of  twelve  honogeneous  linear 
equations  In  c£°^(a  ■  n,  n  +  1),  i.e., 

D  C  -  0,  (19) 

where  D(6)  Is  a  12  x  12  coefficient  natrlx  Involving  &  In  a  transcendental 
form,  and  C  Is  an  unknown  12  x  1  column  eigenvector.  The  nontrivial  solution 
for  C  requires  that  the  determinant  of  the  coefficient  matrix  vanishes,  l.e., 

1D(6)I  -  0.  (20) 

This  leads  to  a  standard  eigenvalue  problem,  and  the  6  can  be  determined  from 
the  transcendental  characteristic  equation.  Standard  numerical  methods  such 
as  the  MUller  method  [22]  with  the  aid  of  a  digital  computer  are  needed  for 
this  purpose.  The  eigenvalues  determined  from  Eq.  20  provide  important 
Information  on  the  fundamental  structure  of  stress  and  displacement  solutions 
for  the  composite  delaminstion  problem.  Furthermore,  the  eigenvalues  6n  which 
satisfy  the  following  condition: 

-1  <  Re [6q]  <  0  (21) 

characterize  the  fundamental  nature  of  stress  singularities  and  provide  the 
asymptotic  stress  and  deformation  fields  at  the  delamination  tip.  In  the  case 
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that  local  crack  surface  tractions  are  nonvanishing ,  for  example,  the  crack 
closure  problem,  Eq.  19  needs  to  be  modified*  Delaminations  having  crack-tip 

•  deformation  configurations  with  nonvanishing  local  traction  boundary 
conditions  are  discussed  in  detail  in  the  next  section. 

*  For  a  delamination  problem  in  composite  laminates  with  general  lamination 

variables  and  fiber  orientations,  the  algebraic  multiplicity  of  the 

eigenvalues  determined  from  Eq.  20  may  give  rise  to  additional  terms  of  the 
6n  lm 

logarithmic  form  Z^  (in  Z^)  in  the  homogeneous  solution,  as  first  suggested 
by  Dempsey  and  Sinclair  [23,24].  In  this  situation,  the  following  terms  may 
also  be  a  part  of  the  homogeneous  solution  in  addition  to  Eqs.  IS  and  16: 

v  6  *m  6 

,  l  Hlk  Z^"]  (i  -  1, 2,4,5, 6),  (22a) 

k"1  36  “ 
n 

»  6  A*m  (6+1) 

“  ’  '  J,  VK  Vs)k  \  (22b) 

lt“1  36  “ 
n 

where  ^  is  the  order  of  the  logarithmic  multiplier  in  the  eigenfunction 
corresponding  to  eigenvalue  6n,  and  is  related  to  Che  property  of  the  D 
matrix  by  ■  M-(N-R),  in  which  M  is  the  algebraic  multiplicity  of  the  root 
6n,  and  N  and  R  are  the  order  and  the  rank  of  the  D  matrix,  respectively.  The 
presence  of  the  logarithmic  terms,  Eqs.  22(a)  and  22(b),  in  the  homogeneous 
solution  equlres  a  nontrivial  solution  for  C^.  Detailed  discussion  of  the 
conditions  for  the  existence  of  Eqs.  22(a)  and  22(b)  in  the  composite 
delamination  problem  can  be  found  in  [25]. 

In  che  construction  of  asymptotic  solutions  for  delamination  stress  and 
displacement  fields,  the  particular  solution  for  the  system  of  governing 
differential  equations  also  contributes  to  the  complete  solution.  It  is 
apparent  that  the  structure  of  the  particular  solution  is  related  to  the 
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u> 


ipplltd  loading  and  detonation  of  the  delaainated  coapoaite.  For  the 
convenience  of  further  developments,  we  consider  here  the  case  of  a  coapoaite 
laalnate  with  delaalnatlons  subjected  to  a  uniform  axial  stretching,  i.e., 
er  ■  eQ.  Under  this  circumstance,  It  has  been  shown  [25]  that  the  particular 

solution  has  a  slallar  fora  as  Eqa.  22(a)  and  22(b), 

•{*’  •  Col  +  <P>  “ll  -0  <l  -  <«■> 

96  »  k-‘ 
n 

“jP>  ■  X  CkP>  "<J+5>k  \°+1/<S«W»)6  -0  <1  ■  (23k> 

d6  °  k_1  11 
n 

where  the  coaponents  ooi  and  uQj  are  known  quantities  determined  by  the  reaote 

loading  condition.  The  *o  in  Eqa.  23(a)  and  23(b)  Is  the  order  of  the 

logarlthaic  eigenfunction  at  5n  -  0  and  la  related  to  the  multiplicity  ^  of 

the  root  6„  »  0  and  the  rank  and  the  order  of  the  aatrix  D  by  !0  *  ^-(N-R) . 

(Note  that  Eqa.  23(a)  and  23(b)  contain  logarlthaic  terms  of  the  forma 

1q 

(In  Zk),  (in  Zk)2,  ...  (in  Zk)  .]  The  necessary  and  sufficient  conditions 
for  the  existence  of  the  particular  solution,  Eqa.  23(a)  and  23(b),  can  be 
shown  [25]  as 

cJ(L)  •  p*  -  0  (24) 

for  every  left  eigenvector  of  D*(0)  defined  in  [25],  where  P*  la  a 

loading  vector  resulting  froa  aol  and  u0j,  and  the  dot  (•)  denotes  the  Inner 

product  of  the  two  coluan  vectors.  In  the  case  that  Eq.  24  does  not  hold,  one 

needs  to  consider  the  logarlthaic  terms  of  a  higher  order  through  a  higher- 

(i +0  (1+1) 

order  differentiation  d  /dfiQ  in  Eqs.  23(a)  and  23(b).  A  detailed 
discussion  on  this  la  given  In  Reference  25. 
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4.  muMmaim  stubs  sukolautiks  n  comfowib  umduiks 

Based  on  cha  general  solution  structures  given  In  Eqs.  15,  16,  22  and  23, 
It  la  possible  at  this  point  to  examine  the  detailed  nature  of  stress 
singularity  associated  with  a  delamination  In  a  fiber-reinforced  composite 
laminate.  Because  of  Che  local  nature  of  Che  stress  singularity,  we  focus  our 
attention  on  the  crack-tip  region  of  a  delamlnation  between  the  mth  and 
(m+l)th  laminae  with  fiber  orientations  9a  and  respectively.  Both  fully 

open  and  partially  closed  delaminatlona  are  considered.  In  the  case  of  a 
delamlnation  with  an  extremely  small  area  of  crack-tip  closure,  a  simplified 
model  by  taking  the  limiting  case  of  a  partially  closed  crack  Is  Introduced. 

4*1  Delamlnation  with  Traction-Free  (Fully  Open)  Surfaces 

Assuming  that  che  crack  surfaces  are  fully  open  and  that  the  interface 
dBj  between  the  piles  Is  perfectly  bonded  along  r  >  0  as  shown  in  Fig.  2,  we 
can  1 Mediately  introduce  the  local  traction-free  boundary  conditions  along 
the  delamlnation  surfaces  *  ■  **, 

o^r.n)  -  ^(r.x)  -  ^(r.x)  -  0,  (25a) 

"  *£fl)<r***>  “  0.  (25b) 

The  continuity  (or  matching)  conditions  of  Interlaminar  stresses  and 
displacements  along  the  ply  interface  6-0  read  as  follows: 

{o<#J,(r,0),T^)(r,0),t^)(r,0)>  -  (o^^r.OJ.t^^^.O.-t^^r.O)),  (25c) 

{u^)(r,0),uJ")(r,0),uj[B)(r,0)>  -  {uJ“+I)(r,0),u^I)(r,0),u^B+l>(r,0».  (25d) 

More  explicitly,  the  local  homogeneous  boundary  conditions,  Eqs.  25(a)  and 
25(b),  and  continuity  conditions,  Eqs.  25(c)  and  25(d),  have  Che  forms  as 


(♦ 
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C7  *e 

of  ?o 


3 

l  (c 

k-1 


<■> 

k 


Hik)(,,)  i4">(1l),fi  +  4+3  Hu)<,‘> 


0. 


(26a) 


H<ri)c-)  + 4+3° 


,(■+1) 

lk 


(-«) 


01®) 


0,  (26b) 


jl««S0«$?  ♦  C<;>  r<;>)  -  ♦  <&*>  T<f  *>..  -  o. 

(i  -  1,2,3;  j  -  1,2, 3, 4, 5, 6). 


where 


4a)<*)  -  +  4a)  •‘u)/o  +  4a>>, 


(26c) 


(26d) 


and 


r<a>-  i  r(a)-  n5a^  r(a)-  u5a^  r(®)«  p(®)«  p(a)« 

llk  l7k  \  ’  l3k  \  '  l4k  pk  »  *5k  qk  ’  r6k  ck 


(a  -  a,  a+1). 


(26e) 


Equations  26(a),  26(b)  and  26(c)  consist  of  a  syater  of  twelve 
hoaogeneous  linear  algebraic  equations  In  4™^  an<*  4"+l>.  The  existence  of  a 
nontrivial  solution  requires  chat  Eq.  20  holds,  leading  to  a  standard 
eigenvalue  problem.  The  solution  for  6n  can  be  obtained  easily  froa  Eq.  20 
and  shown  to  have  the  fora, 

®0  ■  (n  -V2)  *  lY.  (n  -V2),  a-d  n,  (27) 

where  n  -  0,1,2,...*;  and  Y  1»  '  constant  related  to  aaterlal  elastic 
properties,  and  S^*^,  °*  the  adjacent  ath  and  (a+l)th  flber-coaposlte 

laminae.  In  general,  the  value  of  y  needs  to  be  determined  numerically  from 
Eq.  20,  which  Involves  5n  In  a  transcendental  fora.  It  Is  laportant  to  note 
that  Che  eigenvalues  of  An  obtained  froa  Eq.  27  give  critically  laportant 
Information  on  the  fundamental  structure  of  stress  and  displacement  solutions 


for  delaminated  composite  materials.  Wo  mark  Chat  the  eigenvalue* 

(n- V2 )  *  lY  and  (n-  V2  )  are  alngle  root*  and  that  all  the  Integers,  n. 

Including  aero,  have  an  algebraic  multiplicity  of  3  In  general.  Aa  Mentioned 

In  the  preceding  section,  the  velue*  of  6n  which  neet  the  condition  Eq.  21 

provide  the  exact  strength  (or  order)  of  the  Inherently  stress  singularity  for 

the  asynptotlc  stress  solution  at  the  delaal nation  crack  tip.  The  possible 

presence  of  weaker  singularities  and  related  terns  in  logarithmic  fores  as 

1  6  1 

discussed  In  Section  3.3.  i.e.,  (In  Z^)  0  and/or  Z^n(In  Z^)  ",  in  the 
hoaogeneous  solution  as  well  as  in  the  particular  solution  will  be  discussed 
later. 

For  the  purpose  of  llluetration,  consider  a  de lamination  located  along 
the  ply  Interface  of  9/90*  graphite-epoxy  composite  laminae  (Fig.  2).  The 
following  neterlal  elastic  constants*  of  hlgh-aodulus  unidirectional  graphite- 
epoxy  are  used  In  the  computation: 

EL  -  20  x  10*  psi  (137.9  GPa), 

-  Ez  -  2.1  x  10*  pel  (14.48  GPa). 

‘T.T  "  cTz  “  ®Lx  "  0*85  x  10*  psi  (5.86  GPa),  (28) 

VLT  *  vLz  *  *Tz  " 


*Tbese  ply  elastic  constants  are  used  In  the  computation  here  only  to 
Illustrate  the  general  nature  of  the  current  problea.  (These  constants  are 
selected  for  historical  reasons  because  they  have  been  used  In  asny  previous 
studies  of  the  mechanics  of  coaposlte  laminates  [19,26,27].)  Numerical 
results  bssed  on  reel  neterlal  constants  of  the  commonly  employed  T300/5208 
graphite-epoxy  with 

El  -  19.5  x  10*  psi  (134.45  GPa),  Ej  -  Ej  -  1.48  x  10*  pal  (10.2  GPa), 

®LT  “  ®Lx  *  x  10&  P*1  (3.52  GPa),  G-j*  -  0.49  x  10*  psi  (3.38  GPa), 

vtT  "  VL*  *  vTx  *  (28m) 

are  also  given  In  Tables  1  and  2  for  comparison.  The  differences  between  the 
two  cases  are  generally  very  small. 
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where  the  subscripts  L,  T  end  z  denote  the  fiber,  treneveree  end  thlckneee 
directions  of  the  eoaposlta  laaina,  respectively.  The  first  three  elgenveluee 

m 

6t(i  -  1,2,3)  which  satisfy  the  eforeaencloned  constraint  condition  Eq.  21  ere 
given  in  Teble  1  to  illustrete  the  exect  strength  of  the  stress  slngulerity 
essocieted  with  the  delealnetion.  To  demonstrate  further  the  generel 
cherecterlstlcs  of  the  stress  slngulsrltles  for  delealnetion,  results  for  en 
lnterlealner  creek  between  3O*/0  graphite-epoxy  coaposltes  with  the  seae  ply 
properties  ere  elso  shown  (Teble  2)  for  verloue  fiber  orlentetlons  0's.  Froa 
Teblee  1  end  2,  we  observe  thst  e  fully  opened  delealnetion  between  diselailer 
highly  anisotropic  leal nee  always  possesses  three  distinct  stress 
singularities,  l.e.,  e  pair  of  coaplex  conjugates,  6^  2  ■  "  V2*  iy  •  end  e 
reel  constant,  63  “  -0.5.  This  situation  is  unique  end  apparently  different 
froa  the  cases  of  an  interface  crack  between  two  diselailer  isotropic  asdle  or 
orthotropic  solids  in  that  the  three  distinct  dominant  stress  singularities, 
6j,  &2  end  63,  always  exist  slaultaneously  in  the  present  flber-coaposlte 
delealnetion  problea.  In  the  special  cases  when  a  delamination  is  located  in 
the  90*/90*  or  30*/ 30*  coapoaite  system  the  classical  Inverse  square-root 
singularity  for  crack-tip  stresses  is  fully  recovered  as  shown  in  the  Tables, 
because  the  composite  laminate  becomes  unidirectional.  We  note  hare  that  the 
laaglnary  part  of  &i  and  62,  l.e.,  the  value  of  y,  1*  generally  very  saall  as 
compared  with  the  real  part  of  6^  in  all  cases  studied. 

4.2  Pel  sari  nation  with  Ckack-Up  Closure 

,  Froa  Eq.  27  and  froa  the  results  shown  in  Tables  1  and  2,  it  is  clearly 

seen  that  the  asymptotic  delaaination  stress  field  in  dissimilar  anisotropic 
coaposltes  possesses  the  well-known  oscillatory  singularities.  The  associated 
dlsplacesant  field  also  exhibits  oscillatory  characteristics,  leading  to 
controversial  crack-surface  interpenetration  or  overlapping,  which  is 
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physically  inadmissible.  Similar  results  have  also  boon  nocsd  by  several 
lnvestlgstors  In  studying  an  Interface  crack  between  dissimilar  Isotropic 
materials.  In  recent  studies.  Wang  and  Choi  (10.16)  have  shown  that  for  * 
de lamination  between  dissimilar,  strongly  anisotropic  fiber  composites  with 
certain  combinations  of  laminar  elastic  properties,  ply  orientations,  and 
loading  conditions,  global  crack  surface  closure  may  occur.  Under  these 
circumstances,  interlaminar  crack-surface  contact  or  closure  needs  to  be 
considered. 

Consider  the  ease  that  a  da lamination  a  is  located  between  the  ath  and 
(p+l)th  laminae  and  a  portion  of  the  crack  surface,  o,  is  closed  as  shown  in 
Fig.  3.  Frictional  coefficients  associated  with  and  rr^  oo  the  interface 
6  »  0  are  denoted  by  and  /r^,  respectively.  An  exact  analytl&U.  complete 
elaetlclty  solution  for  the  delaai nation  problem  with  crack  closure  is 
generally  difficult  to  obtain  because  the  unknown  contact  stress  distributions 
along  the  crack-closure  region  need  to  be  determined  as  a  part  of  the  final 
solution  (10,14).  However,  the  local  stress  singularities,  asymptotic  field 
solutions,  and  associated  characteristics  can  still  be  detetalned  exactly  by 
using  the  seas  procedure  discussed  in  Section  4.1  but  with  sons 
modifications.  Referring  to  Fig.  3  for  a  partially  closed  delaalnatloo, 
instead  of  using  8qs.  23(a)  and  23(b)  we  can  Introduce  the  local  boundary 
conditions  in  the  crack-closure  region  (-  c  <  r  «*  0)  as  follows: 

^(r.O  -  o^l>(r,-w)  -  0,  T^(r,».  •  ^°(r,-v)  -  0, 

Tr‘>(r'*)  "  ■  °»  *  ■J"*'1>(r.-s)  -  0, 
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T^tr.O  --/M  •|j)(r,n).  (29) 

(a  •  u,  nd). 

Along  thn  ply  Interface  ♦  ■  0,  thn  Man  continuity  conditions  Eqe.  25(e)  nod 
25(d)  (or  Oj°*  and  u*°*  (a  •  a,  aM)  nrn  applicable. 

Using  Eqn.  25  nod  29  nod  following  thn  naan  procedure  given  in  Sect  too 
4.1,  we  can  laaedlately  determine  the  elgenvaluee  4B  for  a  coapoelte 
delaal nation  with  crack-tip  closure.  The  numerical  eiaaple  of  a  delaalnatloo 
located  along  the  Interface  of  9/-9  graphite-epoxy  coepaaltea  la  studied  here 
first.  Stress  singularities  associated  with  the  partially  dosed  delaalnatloo 
crack  tip  with  different  values  of  frictional  coefficients  fT ^  and  are 
shown  In  Table  3.  The  crack-tip  stress  singularity  la  found  to  he  dwnys  -0.5 
with  an  algebraic  eultlplldty  of  2  (l.e.,  double  roots  4j  •  4j  •  -0.5)  (or 
the  delaalnatloo  having  crack  surfaces  in  frictlonlesa  contact  (l.e.,  fr^  • 

0).  In  fact,  the  Inverse  square-root  stress  singularity,  4|  •  4j  * 

-0.5,  la  found  for  all  delaminated  9/-9  fiber  conpoeltea  with  frictlonlesa 
crack-surf  sea  contact.  In  the  cases  of  ft ^  *  0  and/or  /♦.  •  o.  atreaa 
singularities  always  possess  an  Invariant  cooetaat  4|  •  -0.5  (single  root)  as 
In  the  aforementioned  frictlonlesa  contact  cane,  and  a  42  (with  4j  *  -0.5, 
single  root),  dilch  depends  on  values  of  ths  frictional  coeffldenta  (Tables  3 
and  4).  In  Table  3,  values  of  42  for  all  delaalnated  I/-4  fiber  conpoeltea 
studied  are  observed  to  be  slightly  larger  than  -0.5,  wtan  /„  >  0  U 
considered.  That  la,  frictional  contributions  lead  to  a  weaker  delaalnatloo 
stress  singularity  42  than  that  In  a  frictlonlesa  contact  ease  and  In  a 
conventional  homogeneous  open  crack  case.  Ua  note  hern  that  for  a 
delaal net ion  between  9/-9  graphite-epoxy  composites,  the  atreaa  singularities 
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tn  tliwyi  lsdtpMdMt  of  Clio  value  of  the  frlccloool  coofftelooc  fr^ ,  bocouoo 
of  cho  araoetry  (and  antlayuoetry )  of  coopooeata  to  olooctc  aclffneaa  oacrlcaa 
of  Cho  •  aod  -■  pi loo  and  cha  decoupling  of  froo  Oj  and  tyi  la  cho 
forauladoo.  thle  phoooooaoo  la  clearly  aoaa  la  tabla  3,  where  6j  la  always 
HI.)  aad  i]  dlffora  f.»  -0.3  gradually  aa  cho  value  of  tacraaaaa.  NoCe 
fax  char  choc  davladooa  of  <2  ln  the  frlccloaal  contact  caaaa  froo  cha 
ceovendonal  equare-rooc  alagularlcy  ara  rachar  aaall  for  all  6/-8  graphite- 
epoxy  dalxilnacloo  problooo. 

Scroaa  alogularlclaa  ara  alao  determined  for  dalaalnacloos  la  oora 
gaooral  caaaa  of  *|/*2  *lher  coopoolcoo  with  8|  a  82*  r°r  Hiuacraclon,  cha 
raoolca  of  a  dolaaloatloo  bacuooo  30*/90*  graphite-epoxy  eoapoalcaa  ara 
pcaooocad  la  Table  4  for  varloue  valued  of  ff ^  aad  Ic  la  aoaa  froo  cha 

Tahlo  choc  cho  double  root*  6j  •  <2  *  "0*3  alao  appear  for  a  da lamina don  In 
f|/*2  coopoolcoo  wlch  crack  ourfaceo  la  frlcclooloaa  coocacc  (/r^  -  •  0). 

■Mover,  la  cha  caees  of  a  dolaolnaclon  with  crack  eurfacee  la  frlcclonal 
coocacc,  <2  **  dppaconcly  lofluaocod  by  cha  valuoo  of  boch  fr ^  and  In 

table  4,  valuoo  of  <2  *or  the  partially  eloeed  dolaolnaclon  in  30*/90* 
graphlca-apoxy  eoapoalcaa  with  differed  fr^  aad  f^t  ara  obaarvad  to  ba 
aoallar  choa  cha  claaalcal  equare-rooc  acraaa  alagularlcy.  Thua,  62  can  ^ 
etcher  graacar  or  aoallar  than  cho  conventional  lovaraa  equare-rooc 
alagularlcy,  do pending  upon  cha  valuoo  of  aad  f 41  aod  fiber  orlencadona 
at  cho  coopoolcoo.  feting  to  cha  eooplax  algebraic  atrucCura  of  cha 
craoecendencal  choractarlaclc  equation,  lq.  20,  lc  la  generally  00c  poaalbla 
Co  pradlcC  In  explicit  fora  who Char  ^  >  *0.5  or  62  *  "0»3  for  a  dolaolnaclon 
vtch  crack  aurfacac  tn  frlccloaal  contact  without  solving  cha  craoecendencal 
aquation  nuoarlcally.  We  raoork  thac  In  cha  cane  of  a  dolaolnaclon  wlch  crack 
aarfacaa  In  frtcttonleea  contact  between  dlaalollar  anisotropic  oedla,  the 


dominant  stress  singularity,  6  *  -0.5,  has  also  been  determined  Independently 
by  using  a  singular  integral  equation  approach  in  [10,16].  Furthermore,  a 
similar  phenomenon  of  stress  singularity  6  >  -0.5  or  6  <  -0.5  has  been 
observed  in  studying  the  interface  crack  between  dissimilar  isotropic  media 
with  crack  surfaces  in  frictional  contact  [28]. 

4.3  De lamination  with  a  Very  Small  area  of  Crack-Tip  Closure 

The  delamlnation  with  open  crack  surfaces  between  dissimilar  fiber 
composites  has  been  shown  mathematically  in  Section  4.1  to  possess 
controversial  oscillatory  crack-tip  stress  and  displacement  fields.  This 
abnormality  is  thought  to  be  artifacts  resulting  from  the  method  of  approach 
by  using  eigenfunction  expansion  in  the  formulation  and  solution.  As  first 
pointed  out  by  England  [7],  Malyshev  et  al.  [6],  and  later  by  Wang  and  Choi 
[10],  the  region  of  oscillatory  solutions  for  a  delamlnation  with  open  crack 
surfaces  in  a  nominal  tensile  field  is  generally  extremely  small  in  comparison 
with  the  size  of  the  Interface  crack  and  this  very  localized  abnormality  may 
not  be  significant  in  practical  terms  of  linear  fracture  mechanics.  In  fact, 
using  the  partially  closed  crack  model,  Wang  and  Choi  [10]  have  shown  that  a 
composite  delamination  in  a  tensile  field  has  an  extremely  small  crack-tip 
closure  with  c/a  ~  0(10“^).  A  simplified  model  which  disregards  the  small 
closure  (or  oscillatory)  region  and  approximates  the  asymptotic  field  by  an 
inverse  square-root  stress  singularity  is,  therefore,  proposed  for  this 
situation  and  shown  to  provide  excellent  results  comparing  with  those 
determined  by  using  a  partially  closed  crack  model  [10]. 

Under  certain  loading  conditions  other  than  pure  tension,  however,  an 
lnterl  minar  crack  may  also  possess  a  very  small  area  of  crack-tip  closure, 
depending  upon  loading  modes  and  material  elastic  properties  of  the  dissimilar 
composite  laminae  [  16 ] .  A  simplified  solution  for  this  case  can  be  obtained 
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by  taking  the  limit  of  the  crack-cloaure  length,  i.e.,  c  ♦  0,  in  the  reaulta 
derived  from  a  partially  closed  delamination.  In  particular,  the  delamination 
stress  singularity  can  be  taken  directly  from  the  partially  closed  crack  model 
(with  frictionless  crack  surface  contact)  a&  an  inverse  square-root  one. 
Mathematically,  this  is  equivalent  to  finding  an  analytical  solution  for  a 
fully  open  crack  by  following  the  same  formulation  and  procedure  for  a 
partially  closed  crack  case  with  Infinitesimal  closure  length,  and  the 
approximation  Introduced  has  the  effect  of  smoothing  the  oscillatory 
singularity  to  an  Inverse  square-root  singularity  for  the  composite 
delamination.  Therefore,  Interlaminar  stress  intensity  factors  and  strain 
energy  release  rates  can  be  defined  in  a  manner  consistent  with  those  for  a 
homogeneous  crack  and  for  the  refined  model  of  an  interface  crack  between 
dissimilar  isotropic  solids  introduced  by  Comninou  [14].  As  will  be  shown 
later  [18],  this  simplification  leads  to  a  very  effective  and  efficient 
approach  to  the  complex  problem  of  delamination  with  a  fully  open  crack  tip  or 
with  a  very  small  area  of  crack-tip  closure,  and  provides  meaningful 
information  on  the  fundamental  mechanics  of  delamination  problems  in  composite 
laminates  under  general  loading  conditions. 

4.4  Logarithmic  Stress  Singularities 

As  mentioned  in  Section  3  that  besides  the  power-type  stress 
singularities  given  in  Sections  4.1,  4.2  and  4.3  for  various  conditions  in  the 
crack-tip  region,  weak  logarithmic-type  singularities  may  also  appear  in  the 
homogeneous  and  particular  solutions  for  the  delamination  stress  field.  Since 
the  particular  solution  is  related  to  the  remote  loading  applied  to  the 
delaminated  composite,  it  has  to  be  considered  and  constructed  for  each 
individual  case.  To  study  the  possible  presence  of  the  logarithmic  stress 
singularities  in  a  delamination  mechanics  problem,  we  consider  a  symmetric 
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composite  laminate  subjected  to  uniform  in-plane  stretching  with 
ez  -  eq  for  simplicity  and  without  loss  of  generality.  Also,  we  restrict  our 
attention  at  this  point  to  delaminations  located  in  the  following  three 
composite  systems:  9/-9,  9/0*,  and  9/90*  graphite-epoxy  composites.  Based  on 
the  preceding  theoretical  developments  and  the  conditions  for  the  presence  of 
the  logarithmic  terms  given  in  [25],  we  address  each  individual  case 
separately. 

(1)  Delamination  between  9/-9  composites  (9*0*  and  90*) 

Numerical  calculations  by  using  the  ply  elastic  constants  given  in 
Eq.  28  provide  the  following: 

N  -  12,  R  -  9,  M  -  3  (6n  -  integer),  (30a) 

N  -  12,  R  -  11,  M  -  1  <&n  *  integer).  (30b) 

Applying  both  Eq.  30(a)  and  Eq.  30(b)  to  the  condition  1  -  M-(N-R),  we  obtain 

&n  l 

«  0.  Thus,  logarithmic  terms  of  the  form  2^  (*n  Z^)  with  Am  >  1  do  not 

appear  in  the  homogeneous  solution  for  this  class  of  problems. 

Also,  carrying  out  the  computations  of  constructing  the  left 
*(L)  * 

eigenvector  C  and  the  loading  vector  p  ,  we  find  that,  for  all  three  sets 

of  C.  ^  at  &  “  0  In  this  case.  Equation  24  is  satisfied  Identically. 
n  (*+1) 

Therefore,  logarithmic  terms  of  the  form  (in  Z^)  do  not  occur  In  the 

particular  solution  either. 

(2)  Delamination  between  9/0°  or  9/90*  composites  (9  *  0*  and  90°) 

In  these  two  cases,  following  the  same  procedure  and  computations  as 
discussed  in  (I)  but  with  minor  modifications,  we  obtain  similar  results  as 
those  in  the  9/-9  case,  i.e., 

N  ■  12,  R  »  9,  and  M  -  3,  (6n  ■  integer)  (31a) 

N  »  12,  R  -  11,  and  M  -  1,  (5„  *  integer)  (31b) 
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and.  also,  *  0.  Moreover,  Eq.  24  also  holds  for  this  problem.  Thus,  we 

conclude  that  no  logarithmic  singularities  of  any  kind  would  appear  in  the 

asymptotic  solutions  for  delaalnatlon  in  8/90*  and  8/0*  composites;  only 

5n 

power-type  singularities  2^  occur  in  these  problems. 

We  further  remark  that,  In  fact,  it  has  been  shown  in  [25]  that  no 
logarithmic  terms  of  any  kind  would  occur  In  the  solutions  for  a  general  case 
of  a  delamination  located  between  8}  and  82  fiber  composites  with  6j  and  82 
being  any  arbitrary  fiber  orientations. 
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5.  SQMULY  AID  GOKLOSIONS 

The  Mechanics  of  delaalnatlon  in  fiber  coapoalCe  laminates  has  been 
studied.  Formulation  of  the  problem  la  based  on  Lekhnitskii 'a  complex- 
•  variable  stress  potentials  and  basic  relationships  In  laminate  elasticity 
theory  for  anisotropic  fiber  composites.  The  eigenfunction  expansion  method 
used  In  this  study  appears  to  be  a  suitable  approach  to  determine  delaalnatlon 
stress  singularities  and  fundamental  structures  of  stress  and  deformation 
field  solutions.  Stress  singularities  for  a  delamination  are  found  to  be 
related  to  adjacent  ply  material  properties  and  local  traction  boundary 
conditions.  Numerical  results  for  interlaminar  cracks  in  commonly  used 
graphite-epoxy  composites  with  different  fiber  orientations  and  crack-tip 
conditions  are  shown  to  Illustrate  the  basic  nature  of  stress  singularities 
and  general  solutions  for  the  composite  delaalnatlon  problem.  Based  on  the 
Information  obtained,  the  following  conclusions  may  be  reached: 

(1)  Assuming  the  delaalnatlon  Is  fully  open  and  free  from  surface 
traction,  we  find  that  delaalnatlon  stress  singularities  always  possess  an 
oscillatory  form  by  simultaneous  presence  of  three  distinct  eigenvalues, 

-  V2+1Y.  -  V2  “i Y 1  and  -  V2  .  The  oscillatory  stress  singularities  and  field 
solutions  for  composite  delaalnatlon  are  physically  Inadmissible  because  of 
Interpenetration  of  crack-surface  displacements. 

(2)  For  a  delaalnatlon  with  partially  closed  crack  surfaces  in 
frictionless  contact,  the  present  eigenfunction  expansion  approach  always 
gives  an  eigenvalue  6  •  -  V2  with  an  algebraic  multiplicity  of  two  (i.e., 
double  roots),  Indicating  the  classical  square-root  stress  singularity  is 
recovered  in  the  closed  crack  case. 


(3)  In  the  case  of  a  delaai nation  with  crack  surfaces  la  frictional 
contact,  crack-tip  stress  singularities  depend  not  only  on  aaterlal  elastic 
constants  and  fiber  orientations  of  adjacent  plies  but  also  on  frictional 
coefficients  ft ^  and  f^z  along  the  delaaination  surface. 

(4)  The  crack-surface  friction  nay  lead  to  either  a  stronger  or  weaker 
stress  singularity  than  the  conventional  inverse  square-root  one,  depending 
upon  fiber  orientations  of  the  adjacent  plies.  Present  numerical  results,  for 
exaaple,  show  chat  a  weaker  stress  singularity,  i.e.,  0  >  5  >  -  V2.  occurs  for 
a  delaaination  between  any  8  and  -8  fiber  coaposites,  but  a  stronger 
singularity,  i.e.,  -V2  >  6  >  -1,  occurs  for  a  delaaination  between  30*  and 
90*  coaposites,  if  /r^  >  0  and  f^z  >  0. 

(5)  In  the  situation  that  the  delaaination  contains  a  very  snail  area  of 
crack-surface  closure  (e.g.,  c  <  10  *a),  a  siaplified  model  with  the  crack-tip 
stress  field  having  an  Inverse  square-root  stress  singularity,  as  deternined 
by  finding  the  solution  from  the  Halting  case  of  a  partially  dosed  crack 
solution,  is  suggested  and  later  used  for  solving  the  coaplete  boundary  value 
problea. 

(6)  Examining  the  aultlplicicy  of  eigenvalues  and  the  rank  and  order  of 

Che  coefficient  aatrlx  In  the  eigenfunction  solution,  we  find  that  no 

logarithmic  stress  singularities  of  any  kind  would  appear  in  the  hoaogeneous 

and  particular  solutions  for  the  composite  delaaination  prob1  a;  only  power- 

«n 

type  singularities  of  the  fora  could  occur. 

(7)  After  determining  of  all  the  eigenvalues  for  each  individual 
delaaination  case,  general  solution  structures  for  coaposlte  deformation  and 
stress  fields  can  be  established  laaediately.  Numerical  aethods  such  as  the 
singular  flnite-eleaent  technique,  which  can  incorporate  exact  delaaination 
stress  singularities  in  the  element  formulation,  can  be  easily  developed  to 
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(♦ 


solve  chs  complete  boundary-value  problea  for  delaalnaCions  in  composite 
laal nates  with  any  arbitrary  combinations  of  lamination,  geoaetric,  and  crack 
variables.  One  of  such  aethods  employing  displaceaent-based  singular  crack- 
tip  eleaents  is  given  in  the  associated  paper  [18]. 
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TABLE  1 

DOMINANT  STRESS  SINGULARITIES  *  FOR  DELAMINATION 
BETWEEN  9/90*  GRAPHITE-EPOXY  COMPOSITES 


e 

3I 

*2 

*3 

o*+ 

-0.5 

+ 

o.osnoi: 

-0.5 

0. 05110 1-: 

-0.5 

(-0.5 

+ 

0.032924-:) 

(-0.5 

- 

0.032924-:) 

(-0.5) 

15* 

-0.5 

+ 

0.050349*: 

-0.5 

• 

0 .050349t 

-0.5 

(-0.5 

+• 

0.032484*) 

(-0.5 

- 

0.032484- > 

(-0.5) 

30* 

—0 .  o 

+ 

0.045138*: 

-0.5 

0.045138t 

-0.5 

(-0.5 

+ 

0.025764*) 

(-0.5 

- 

0.025764*:) 

(-0.5) 

45’ 

-0.5 

♦ 

0. 034504  •: 

-0.5 

• 

0.034  504-: 

-0.5 

(-0.5 

+• 

0.015604*:) 

(-0.5 

- 

0.015604*) 

f-0.5) 

60* 

-0.5 

4- 

0.021U9-: 

-0.5 

_ 

0.021119*: 

-0.5 

(-0.5 

4- 

0.008067*) 

(-0.5 

- 

0.008C67*:) 

(-0.5) 

75* 

-0.3 

+ 

0.008899*: 

-0.5 

• 

0.008899: 

-0.5 

(-0.5 

4- 

0.004265*:) 

(-0.5 

- 

0.004265*:) 

(-0.5) 

90’" 

-0.5 

-0.5 

-0.5 

r-0.5) 

(-0.5) 

(-0.5) 

*Yalues  la  the  parentheses  are  for  1300.’  5280  graph  Ice-epoxy  with  Laminar 
elastic  constants  given  In  Eq.  28(a). 

■A 

These  eigenvalues  are  for  0*/90*  and  90’/90*  composites  in  a  general  loading 
condition.  In  the  cases  of  0’/90’  and  90'/90*  composites  under  uniform 
stretching  e_  ■  sQ,  63  »  -0.5  does  not  appear  because  of  :u  ■  r  •  0  and 
being  decoupled  from  ocher  stress  components.  y 
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TABU  2 

DOMINANT  STUBS  SINCOUBITIES*  FOB  DELAMI NATION 
BETWEEN  30*/«  OUreiTT -EPOXY  COMPOSITES 


e 

*1 

*2 

*3 

0* 

-0.5  ♦ 

0.0124511 

-0.5  - 

0.012451i 

-0.5 

(-0.5  ♦ 

0. 009456t) 

(-0.5  - 

0.0094560 

(-0.5) 

15* 

-0.5  ♦ 

0.0 1049 It 

-0.5  - 

0. 01049U 

-O.S 

(-0.5  ♦ 

0.0079200 

(-0.5  - 

0.0079200 

(-0.5) 

30* 

-0.5 

-0.5 

-0.5 

(-0.5) 

(-0.5) 

(-0.5) 

45* 

-0.5  ♦ 

0.01596M 

-0.5  - 

0.0159681 

-0.5 

(-0.5  ♦ 

0. 011608i) 

(-0.5  - 

0.0116080 

(-0.5) 

60* 

-0.5  ♦ 

0.030943t 

-0.5  - 

0.0309431 

-0.5 

(-0.5  ♦ 

0. 020441i) 

(-0.5  - 

0.0204410 

(-0.5) 

73* 

-0.5  ♦ 

0.041030t 

-0.5  - 

0.0410301 

-0.5 

(-0.3  + 

0.0246640 

(-0.5  - 

0.0246640 

(-0.5) 

90* 

-0.5  ♦ 

0.045138* 

-0.5  - 

0.043138i 

-0.5 

(-0.5  ♦ 

0.0247640 

(-0.3  - 

0. 025764i) 

(-0.5) 

*Valuas  la  tht  parent haaas  an  tor  T300/5280  graph lt*-«poxy  with  Ualatr 
•laatlc  propart glvaa  la  Eg.  28(a). 
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TABLE  3 

DOMINANT  STUBS  SINCULAUTIES*  POK  DC LAKI NATION 
WITH  C1AOC-TIP  CLOSUtt  IN  •/-•  OUPHITE -EPOXY  COMPOSITES 


15* 

30* 

45* 

60* 

75* 

0.0 

-0.5 

-0.5 

-0.5 

-0.5 

-0.5 

-0.5 

-0.5 

-0.5 

-0.5 

-0.5 

0.2 

-0.5 

-0.4904 

-0.5 

-0.4962 

-0.5 

-0.4940 

-0.5 

-0.4941 

-0.5 

-0.4964 

0.4 

-0.5 

-0.4977 

-0.5 

-0.4923 

-0.5 

-0.4661 

-0.5 

-0.4681 

-0.5 

-0.4928 

0.6 

-0.5 

-0.4965 

-0.5 

-0.4664 

-0.5 

-0.4621 

-0.5 

-0.4822 

-0.5 

-0.4892 

o.a 

-0.5 

-0.4954 

-0.5 

-0.4646 

-0.5 

-0.4762 

-0.5 

-0.4763 

-0.5 

-0.4855 

*«!  and 

<2  found 

co  b«  lndopmdant 

of  tho  vo loo  of 

/ 
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TAB  LIC  4 

DOMINANT  STRKSS  SINGULARITIES  FOR  DELAMINATION 
WITH  CRACK-TIP  CLOSURE  IN  30°/90°  GRAPHITE /EPOXY  COMPOSITES 


0.0 

0.1 

0.3 

0.5 

0.7 

-0.5 

-0.5 

-0.5 

-0.5010 

-0.5 

-0.5031 

-0.5 

-0.5051 

-0.5 

-0.5072 

-0.5 

-0.5036 

-0.5 

-0.5046 

-0.5 

-0.5067 

-0.5 

-0.5087 

-0.5 

-0.5108 

-0.5 

-0.5108 

-0.5 

-0.5118 

-0.5 

-0.5138 

-0.5 

-0.5159 

-0.5 

-0.5179 

-0.5 

-0.5179 

-0.5 

-0.5189 

-0.5 

-0.5210 

-0.5 

-0.5230 

-0.5 

-0.5251 

-0.5 

-0.5251 


-0.5 

-0.5261 


-0.5 

-0.5281 


-0.5 

-0.5302 


-0.5 

-0.5322 


7.  LIST  OF  FIGURE  CAPTIONS 


Delaminatlons  in  a  (9^/62/03/. . ./83/82/6X)  Fiber-Reinforced  Composite 
Laminate. 

Coordinates  and  Geometry  of  a  Delamination  with  Open  Crack  Surfaces 
between  9a  and  8^+^  Plies. 

Coordinates  and  Geometry  of  a  Delamination  with  Finite  length  of 
Crack-Surface  Closure  between  9^  and  8^^  Plies. 


*  ‘VW— /VV«ll  FIb«r-Rel„force3  Co.po.lt. 
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OF  POOR  QUA-iVY 

APPENDIX  1 

Expressions  for  coefficients  H^(4)  i°  E<13* 

IS  and  16  are  aa  follows 

Hlk  “  ^uitsln^  +  c°8l>)2* 

H2k"  “nk/pk8ln*  +  co8^ 

■  -(u^siwfr  +  cosAXu^cos*  -  sim>). 

H4k  “  -  ain^)2. 

H5k  *  nk(ufccos4  -  sin*). 

H6k  ”  Pkcos*  +  qk*i»6. 

H7k  "  _Pk8ln^  +  qkcos** 

'’8k  "  ck» 

where  and  are  defined  in  Eqs.  11#  12  and  13,  and  p^,  and  are 
complex  constants  relaced  co  laminar  elastic  constants  by 

Pic  ■  l\\vl  +  2u  "  314\  +  2lsVk  *  2lbV 

\  ■  2l2Pk  +  \zJvV.  "  Vi\/\  +  225\  "  *26’ 

\  *  214pk  +  224/pk  “  244nk/Mk  +  245nk  “  246* 


(A-l) 


(A-2) 


